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Using a sample of 14 million ip(2S) events collected with the BESII detector, branching 
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fractions or upper limits on the branching fractions of ip(2S) decays into 7pp, ~{2(ir + ir ), 
jK^K+tt- +c.c, -yK+K-iv+TT-, <yK*° R-iv + + c.c, jK^'K* , -yn+TT-pp, j2(K + R-), 73(tt + 70 
and j2(tt + -k~)K + K~ with hadron invariant mass less than 2.9 GeV/c 2 are reported. We also report 
branching fractions of ip(2S) decays into yppir , ppTr ir , 2(tt + ty^)ty" , ojtt + tt~ , cj/2(1270), bf 7r T , 
n°K%K + n- +c.c, I&p^IQ, iv°2(-w + n~)K + K~ and 1 -k°2{tt + k-)K + K~ . 

PACS numbers: 13.20.Gd, 12.38.Qk, 14.40.Gx 



I. INTRODUCTION 

Besides the conventional meson and baryon states, 
QCD (Quantum Chromodynamics) also predicts a rich 
spectrum of so-called QCD exotics, among them glue- 
balls (gg), hybrids (qqg), and four-quark states {qqqq) in 
the 1.0 to 2.5 GcV/c 2 mass region. Therefore, the search 
for evidence of these exotic states plays an important role 
in testing QCD. Radiative decays of quarkonium are ex- 
pected to be a good place to look for glueballs, and there 
have been many studies performed in J/tp radiative de- 
cays [1, 2], while such studies have been limited in tp(2S) 
radiative decays due to low statistics in previous experi- 
ments [2, 3]. Radiative decays of tp(2S) to light hadrons 
are expected at about 1% of its total decay width [4]. 
However, the previously measured channels only sum up 
to 0.05% [3]. 

In this paper we present measurements of tp(2S) decays 
into jpp, 72(7r + 7r~), jKgK + ir~ + c.c, ~fK + K~ir + ir~ , 

1 K* a K-Tr+ + c.c, jK* K*°, -fir+TT-pp, j2(K+K~), 
73(7r + 7r~) and -f2(7r + ir~)K + K~ with the invariant mass 
of the hadrons (rrihs) in each final state less than 
2.9 GeV/c 2 . Measurements of ip{2S) decays into 
7PP7T , ppTr°ir , 2(7T+7r-)7r , WTT+Tr", w/ 2 (1270), bfw* 
, n°K%K + ir- + c.c, ^p^K®, Tr 2(ir + ir-)K+K- and 
jTr°2(ir + ir~)K + K~ are also presented and are used for 
the background analysis. Many of the above measure- 
ments have been published previously [5] ; this paper pro- 
vides more detailed information. 



II. BES DETECTOR 

BES is a conventional solenoidal magnetic detector 
that is described in detail in Refs. [6, 7]. A 12-layer ver- 
tex chamber (VTC) surrounding the beam pipe provides 
trigger and track information. A forty-layer main drift 
chamber (MDC), located radially outside the VTC, pro- 
vides trajectory and energy loss (dE/dx) information for 
charged tracks over 85% of the total solid angle. The 
momentum resolution is a p /p = 1.78%-y/l + p 2 (p in 
GeV/c), and the dE/dx resolution for hadron tracks is 
~ 8%. An array of 48 scintillation counters surrounding 
the MDC measures the time-of-flight (TOF) of charged 
tracks with a resolution of ~ 200 ps for hadrons. Radially 
outside the TOF system is a 12 radiation length, lead-gas 
barrel shower counter (BSC). This measures the energies 
of electrons and photons over ~ 80% of the total solid 
angle with an energy resolution of a E /E = 21%/y/E (E 



in GeV). Outside of the solenoidal coil, which provides 
a 0.4 Tesla magnetic field over the tracking volume, is 
an iron flux return that is instrumented with three dou- 
ble layers of counters that identify muons of momentum 
greater than 0.5 GeV/c. 

The data sample used in this analysis was taken with 
the BESII detector at the BEPC storage ring at an 
energy of y/s = 3.686 GeV. The number of ip{2S) 
events is (14.0 ± 0.6) x 10 6 , determined from inclusive 
hadronic events [8], and corresponds to a luminosity of 
£3.686 = (19.72 ±0.86) pb _1 [9], measured with large 
angle Bhabha events. Continuum data, used for back- 
ground studies, were taken at i/s — 3.650 GeV with a 
luminosity of £ 3 .650 = (6.42 ± 0.24) pb _1 [9]. The ratio 
of the two luminosities is £3.686 /£3.650 = 3.07 ± 0.09. 

Monte Carlo (MC) simulations are used for the de- 
termination of mass resolutions and detection efficien- 
cies, as well as background studies. The simulation of 
the BESII detector is GEANT3 based, where the inter- 
actions of particles with the detector material are sim- 
ulated. Reasonable agreement between data and Monte 
Carlo simulation is observed [10] in various channels such 
as e+e~ — > {"{)e + e~ , e+e~ — ► (■y)p + p^, J/tp — > pp and 
i>(2S) -» n+TT-J/ij, J/tp -> £+£- (£ = e,fi). An in- 
clusive ip{2S) decay MC sample of the same size as the 
ip{2S) sample is generated by LUNDCHARM [11] and 
used to estimate backgrounds, 

III. EVENT SELECTION 

A neutral cluster is taken as a photon candidate when 
the following conditions are satisfied: the energy de- 
posited in the BSC is greater than 50 MeV; the first hit is 
in the beginning six radiation lengths; the angle between 
the cluster and the nearest charged track is greater than 
15°; and the difference between the angle of the cluster 
development direction in the BSC and the photon emis- 
sion direction is less than 37°. 

Each charged track is required to be well fitted to 
a three-dimensional helix, be in the polar angle region 
I cos6>| < 0.8 in the MDC, and have a transverse momen- 
tum greater than 70 MeV/c. The particle identification 
chi-squared, Xf/u(*)j is calculated based on the dE/dx 
and TOF measurements with the following definition 

X%i D (i) = XdE/dx(*) + Xtof(*)- 

For all analyzed decay channels, the final states of 
the candidate events must have the correct number of 
charged tracks with net charge zero. 
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If there is more than one photon candidate in 
an event, the candidate with the largest energy de- 
posit in the BSC is taken as the radiative photon 
in the event, and a four-constraint kinematic fit (4C- 
fit) is performed. The combined confidence level, 

prob( X 2 comb ,ndf) = £L t ^U^dt, is 

required to be greater than 1%, where ndf is the num- 
ber of degrees of freedom and x 2 omb ^ s defined as the 
sum of the x 2 01 the kinematic fit (xlc) an d X%id : 
Xlomb = Xlc + E Xp/d(*)> whcrc i runs ov er all charged 

i 

tracks. 

For each decay mode, rrihs is required to be less than 
2.9 GeV/c 2 to exclude ip(2S) radiative transitions into 
other charmonium states, such as XcJi J/ipy an d Vc- To 
remove background from charged particle misidentifica- 
tion, the value of xl om b f° r — > 7 + hs is required 

to be less than those for ip(2S) decays into background 
channels 7 + hs' , where hs' has the same number of 
charged tracks as hs, but different particle types for the 
charged tracks. If there is potential background from 
i>(2S) — ► 7r + 7r~ J/ip, it is largely suppressed by applying 

\ m recZil ~ m JlA > °- 05 GeV/c 2 , where m^Vj is the 
mass recoiling from each possible tt + it~ pair. Possible 
background from ip(2S) — ► ifg + X is removed by requir- 
ing that the invariant mass of 7r + 7r~ is outside the Kg 
mass region (Im^ — m K o | > 0.04 GeV/c 2 ). 



be the normalization factors, but the weights of the sig- 
nal (w s ) and the other backgrounds (wb) are free. With 
this method, the number of signal events is extracted for 
each radiative decay mode for to^ s < 2.9 GeV/c 2 . 
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FIG. 1: The fitted xlc distribution for ip(2S) -> j2(tt + 7T~) 
candidate events. The dots with error bars are data. The 
solid line is the fitted result, which is the sum of the four 
components: signal events (dashed line), MC simulated multi- 
photon backgrounds (dotted line), continuum (hatched his- 
togram), and other backgrounds (dot-dashed line). 



V. EVENT ANALYSIS 



IV. BACKGROUNDS AND FITTING 
PROCEDURE 

In our analyses, the backgrounds for each ip(2S) — > 
7 + hs decay mode fall into three classes: (1) continuum 
background, estimated using the continuum data; (2) 
multi-photon backgrounds, e.g. ip(2S) — > 7r° + hadrons, 
3~f+hadrons, etc., where hadrons have the same charged 
tracks as the signal final state, estimated with MC simu- 
lation and normalized according to their branching frac- 
tions; and (3) other backgrounds, estimated using an in- 
clusive ip(2S) MC sample of 14 million events [11]. Multi- 
photon backgrounds are dominant; continuum back- 
ground and other backgrounds including contamination 
between the channels studied are lower. The observed 
x\c distributions include both signal events and all of 
these backgrounds. 

The number of signal events for most radiative decay 
channels is extracted by fitting the observed x\c distri- 
butions with those of the signal and background chan- 
nels [12], i.e. xlu = WsXlig + T, Wb . w bi xl g , where w s 
and are the weights of the signal and the background 
decays, respectively. As an example, Fig. 1 shows the ob- 
served x 2 distribution for ip(2S) — > j2(ir + Tr~), together 
with the x 2 distributions for the signal, multi-photon, 
continuum, and other background channels, as well as the 
final fit. In the fit, the weights of the multi-photon back- 
grounds and the continuum backgrounds (wb) are fixed to 



A. ip(2S) — ► ~/pp 

The major backgrounds to ip(2S) — > jpp come from 
the channels ip(2S) — > ir°ir°pp, r yTr°pp, and ir°pp. In or- 
der to estimate these backgrounds, we first measure their 
branching fractions. 



1. Background estimation 

For ip(2S) — > ir a iT a pp, candidate events must have two 
charged tracks and four good photons. To reject back- 
ground from ip(2S) — ► 7r°7r° J/ip, J/ip — > pp events, the 
pp invariant mass is required to satisfy \m pp - — mj/^,\ > 
0.1 GeV/c 2 . Figure 2 (a) shows the scatter plot of 
7172 versus 7374 invariant mass for events after selec- 
tion, where 7172 and 7374 arc formed from all possi- 
ble combinations of the four photon candidates. The 
cluster of events shows a clear 7r°7r° pair signal. Fig- 
ure 2 (b) shows the 7172 invariant mass after requir- 
ing the other two photons be consistent with being a 
7T° (|m 7374 — m lx o\ < 0.03 GeV/c 2 ). A fit to the peak is 
performed with a double Gaussian function, where the 
parameters are determined from MC simulation, plus a 
second order polynomial to describe the smooth back- 
ground. The number of events in the peak determined 
from the fit is 254 ± 24, and the background contamina- 
tion to the peak is estimated to be 51 ±11 from the ti 7374 
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FIG. 2: (a) Scatter plot of 7172 versus 7374 invariant mass for 
ip(2S) — ► ty°iy°pp candidate events, (b) Invariant mass distri- 
bution of 7172 after requiring |m 7374 — m^o| < 0.03 GeV/c 2 . 
The data are fitted with a double Gaussian function and a 
second order background polynomial. 



FIG. 3: Distributions of 77 invariant mass for ip(2S) — > 777PP 
candidate events fitted with a double Gaussian function and 
a second order background polynomial, (a) signal and (b) 
backgrounds including continuum, known simulated back- 
grounds (tp(2S) — > TY°iT pp and ir°pp), and other unknown 
backgrounds estimated from the ip(2S) inclusive MC sample. 



sidebands: (0.0, 0.06) and (0.2, 0.26) GeV/c 2 . No signif- 
icant resonance is observed in the mass distributions of 
any possible combination of particles in this final state, 
indicating no significant intermediate processes in the ob- 
served '0(25') — > TT TT°pp candidate events. Therefore the 
detection efficiency for ip{2S) — > 7r°7r°pp is determined to 
be 8.3% using a phase space generator, and the branching 
fraction is determined to be 

B(ip(2S) -> ir°ir°pp) = (1.75 ± 0.21) x 10~ 4 , 

where the error is statistical. 

For ip(2S) —>■ j^pp, candidate events must have two 
charged tracks and three good photons. To reject back- 
ground from ip(2S) — + 7r°7r° J/tp, J/ip — > pp events, the 
pp invariant mass is required to satisfy \m P p — mj/^\ > 
0.1 GeV/c 2 . 

Figure 3 (a) shows the 77 invariant mass distribution 
for ip(2S) — > 777PP candidate events, where 77 is any 
possible combination among the three photon candidates. 
There is a clear 7r° signal. The distribution is fitted with 
a double Gaussian function with parameters determined 
from MC simulation plus a second order polynomial for 
the background. The number of events determined from 
the fit is 345 ± 33. 

Background studies indicate that the main contami- 
nation to the 7T° signal comes from ip(2S) — > 7r°7r°pp 
and TT°pp; other backgrounds only contribute a smooth 
background. All possible backgrounds, including contin- 
uum, known simulated backgrounds (ip(2S) — > ir°Tr°pp 
and Tr°pp), and other unknown backgrounds estimated 
from the ij)(2S) inclusive MC sample, are combined in 
Fig. 3 (b). Fitting this distribution in the same way as 
in Fig. 3 (a), the number of peaking background events 
is estimated to be 219 ± 18. 

Just as for ip(2S) — * ir°ir°pp, no significant intermedi- 
ate process is observed in ?p(2S) — > jir°pp. The efficiency 



is determined to be 8.94% with a phase space generator, 
and the branching fraction is calculated to be 

B(ip(2S) -> jTT°pp) = (1.0 ± 0.3) x 10~ 4 , 

where the error is statistical. 

For %l)(2S) — > tt°pp, there is an earlier measurement 
from BESII [13]. We reanalyze this channel in the same 
way, and then extract the m p p distribution and estimate 
its contamination to i/)(2S) — » 7pp. 

Candidate events are required to have two charged 
tracks and two good photons. The probability of the 4C- 
fit must be greater than 1%, and the probability of the 
4C-fit for the ip(2S) — » 77PP hypothesis must be greater 
than that for ip(2S) — » jjK + K~. To reject background 
from ip(2S) — > 7T 7r° J/tp, J/ip —> PP events, the invariant 
mass of pp is required to be: \m pp — mj/^\ >0.02GeV/c 2 . 

A fit to the distribution is performed with a dou- 
ble Gaussian function with parameters determined from 
MC simulation plus a second order polynomial for the 
background for ip(2S) — > "fjpp candidate events. The 
number of events determined from the fit is 266 ±20, and 
the detection efficiency is 14.8%. The branching fraction 
is determined to be: 

B(ip(2S) -> 7T°pp) = (13.0 ± 1.0) x 10~ 6 , 

where the error is statistical. This measurement agrees 
well with the previous BESII result of (13. 2± 1.0± 1.5) x 
10~ 6 [13]. 

2. Signal Analysis 

For ip(2S) — > 7pp, 329 events are observed after the 
event selection described in Section III; the pp invariant 



■5 



mass distribution is shown in Fig. 4. After subtracting 
the normalized major backgrounds, ?p(2S) — > n n pp, 
pp, and ir°pp, the number of signal events is 142 ± 18. 
The detection efficiency determined from MC simulation 
is 35.3%, and the branching fraction for this process is 
determined to be: 

B(ip(2S) -» jpp) = (2.9 ± 0.4) x 1(T 5 , 

where the error is statistical. 

There is an excess of events between pp threshold and 
2.5 GeV/c 2 , but no significant narrow structure due to 
the X(1859), that was observed in J/tp — > jpp [14]. A 
fit to the mass spectrum (see Fig. 5) with an acceptance- 
weighted S'-wave Breit-Wigner for the X resonance (with 
mass and width fixed to 1859 MeV/c 2 and 30 MeV/c 2 , 
respectively), together with the normalized MC his- 
tograms for the above measured background channels 
(tp(2S) — ^ Tr°Tr°pp, 77T°pp, and Tr°pp) and the histogram 
from ip(2S) — > -fpp phase space [15], yields 11.7 ± 6.7 
events with a statistical significance of 2.0 a. The upper 
limit on the branching fraction is determined to be 

B(ip{2S) -> 7 A(1859) -» 1P p) < 5.4 x 10~ 6 

at the 90% C.L. 
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FIG. 4: The pp invariant mass distribution for tp(2S) — > "/pp 
candidate events (dots with error bars). The shaded his- 
togram is the sum of all backgrounds, including continuum, 
known simulated backgrounds (tp(2S) — ► 7r°7r°pp, 77r°pp, and 
7r°pp), and other unknown backgrounds estimated from the 
ip{2S) inclusive MC sample. 



B. ip(2S) -> 72(77+77-) 

For i[)(2S) — * 72(7r + 7r~), the main background comes 
from tp(2S) — > 2(7r + 7r")7r°, so we first measure tp(2S) — > 
2(ir + TT~)ir° in order to be able to estimate its contami- 
nation to i>(2S) -> 72(7r+7r"). 
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FIG. 5: The fit to the m P p — 2m p distribution of tp(2S) — > 
"fpp candidate events. The solid histogram is the fit result, 
the lower dashed line is the X resonance shape, the dash- 
dotted histogram is the shape for ip(2S) — > "/pp phase-space, 
the dotted histogram is the measured background channels 
(ip(2S) — > n°ir°pp, fw°pp, and TT°pp), and the top dashed line 
is the efficiency curve. 



1. Background analysis 

For tp(2S) — ► 7r°2(7r + 7r _ ), candidate events are re- 
quired to have four charged tracks and two good pho- 
tons. The probability of the 4C-fit must be greater than 
1%, the xLnfc probability for the i/)(2S) -» 7r°2(7r+7r-) 
hypothesis must be greater than those of "4>(2S) — > 
ir°K + K~ir + ir~ and tp(2S) — ► 7r°7r + 7r _ pp, and the sum 
of the momentum of any 7r + and 7r _ pairs must be 
greater than 550 MeV/c to reject contamination from 
ijj(2S) — ► ir + -K~ J/ip events. 

After the above selection, a clear 7r° signal can be 
seen in the 77 invariant mass distribution of ip(2S) — > 
2(7r + 7r _ )77 candidates. After subtracting backgrounds, 
such as ^(25) -> 7X c o,XcO -» 2(7r+7r"), ^(25) -> 
7r°7r° J/V', tZ/V' 2(7r + 7r _ ), etc., in the 77 invariant mass 
spectrum, the distribution is fitted with a 7r° signal shape 
determined with MC simulation plus a second order poly- 
nomial for the other remaining backgrounds, and the 
number of 7r° signal events is 2173 ± 53. The detection 
efficiency is determined to be 6.32% taking into consider- 
ation the significant intermediate states such as oj7r + 7r~ , 
w/2(1270), and b 1 tt t described below. 

Figure 6 (b) shows the 7r + 7r~7r° invariant mass distri- 
bution for events satisfying |m 77 — m n o\ < 0.03 GeV/c 2 . 
The fit is performed with an u> signal shape plus a second 
order polynomial for the background, and the number of 
lo signal events obtained is 386 ±23. The efficiency deter- 
mined from MC simulation is 3.74% correcting for inter- 
mediate states, such as o>/ 2 (1270) and bfir T , described 
below. 

Figure 6 (c) shows the distribution of n + ir~ invariant 
mass recoiling against the u>, selected with the require- 
ments iTOTr+Tr-Tr — < 0.05 GeV/c 2 and |m W7r — raftj > 
0.2 GeV/c 2 to reject bf-rT T events. The invariant mass 
spectrum is fitted with a a, a / 2 (1270) shape determined 
from MC simulation, and a second order polynomial to 
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describe other backgrounds. The number of /2(1270) 
events obtained is 57 ± 13, and the detection efficiency 
determined from MC simulation is 3.65%. 

Figure 6 (d) shows the um* 1 invariant mass spectrum 
with the requirement Im^+^-^o — m u \ < 0.05 GeV/c 2 . 
A clear signal is seen. Fitting with a b\ signal shape 
with the mass and width fixed to PDG values plus a 
background polynomial, the number of bf 1 signal events 
is 202 ± 21, and the detection efficiency determined from 
MC simulation is 3.24%. The branching fractions of these 
processes are determined to be: 

B(ip(2S) -» vr 2(7r+7r-)) = (24.9 ± 0.7) x 10~ 4 , 
B(iJj(2S) -> lutt+tt-) = (8.4 ± 0.5) x 10~ 4 , 
B(ip(2S) -> w/ 2 (1270)) = (2.3 ± 0.5) x 10~ 4 , 
B(^(2S) bf tt t ) = (5.1 ± 0.6) x 10~ 4 , 

where the errors are statistical. 



pion. The background from ip(2S) — > it + it J/ip is re- 
jected by requiring |mJT ec „ iZ — mj/^\ > 0.05 GeV/c 2 . Af- 
ter selection, 1697 candidates remain, and the 2(tt + 7t~) 
invariant mass distribution for the candidate events is 
shown in Fig. 7. The backgrounds include contributions 
from the continuum (estimated from the data sample at 
y/s = 3.65 GeV), tp(2S) -> 7r°2(7r + 7r-), backgrounds re- 
maining from 7r + 7r~ J/ip, J/tp — > pTr and tt KgK + TT~ + 
c.c, and the other unknown backgrounds assuming that 
they have the same shape as that obtained from the in- 
clusive ip(2S) MC sample. 

Using the % 2 fitting method of Section IV, the number 
of signal events is 583±41. The detection efficiency deter- 
mined from MC simulation is 10.4%, and the branching 
fraction for this process is determined to be: 

B(<ip(2S) -> 727T+7T-) = (39.6 ± 2.8) x 10~ 5 , 

where the error is statistical. 
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FIG. 7: The 2(n + n~) invariant mass distribution for 
ip(2S) — ► 72(7r + 7r~) candidate events (dots with error 
bars). The shaded histogram includes contributions from 
the continuum (estimated from the data sample at y/s = 
3.65 GeV), 4>{2S) — ► 7r°2(7r + 7r _ ), backgrounds remaining 
from tt + it~ J/ip, J/tp — > pir and ir°K%K + ir~ + c.c, and other 
unknown backgrounds estimated from the inclusive ip(2S) 
MC sample. 



FIG. 6: Invariant mass distributions with fits for ip(2S) — > 
7r°2(7r + 7r _ ) candidates, where dots with error bars are data, 
and the solid histograms and curves are the fit results, (a) 77; 
(b) 7T + 7i-~7r° with |m 77 -m n o\ < 0.03 GeV/c 2 ; (c) tt + tt~ with 
[m^+^-^a — m w \ < 0.05 GeV/c 2 and bfn^ events rejected; 
and (d) ojtt* for the ip{2S) -> 7r°2(7r + 7r-) candidate events. 
Resonance parameters are fixed to their world averaged val- 
ues [3]. 



C. ip(2S) -> ~iK%K+tt- + c.c. 

For ip(2S) — > jKgK^TT^ , the main background comes 
from ip(2S) — > K s K ± tt t tt° , so we first measure ip(2S) — > 
KgK^n^n in order to estimate its contamination to 



2. Signal analysis 

For ip(2S) — > 72(7r + 7r~), candidate events require four 
charged tracks, and each track must be identified as a 



1. Background estimation 

For ip(2S) — > 7r°i ; ir^i ; sr + 7r _ + c.c, candidate events re- 
quire four charged tracks and two good photons. Fig- 
ure 8 shows the scatter plot of 7r + 7r~ invariant mass ver- 
sus the decay length in the transverse plane (L xy ) of Kg 
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candidates, where a clear K s signal is observed. Can- 
didate events are required to have only one K s can- 
didate satisfying the requirements \m w + v - — m K o\ < 
0.015 GeV/c 2 and L xy > 0.5 cm. After K° s selec- 
tion, the remaining two tracks are identified using their 

the final state 



X w +k- 



if x\. 



< 



K+> 



where 



Xk-k values, i.e., if X K + E ~ < 
is considered to be j^K s K + 7r 
the final state is considered to be jK s K~n~ 

= Xpid( K ) + XpidW- 

The confidence level of the 4C-fit must be greater than 
1%, and the sum of the momentum of any 7r + and 7r~ pair 
greater than 550 MeV/c to reject contamination from 
ip(2S) — ► Tt + iT~ J/ip events. 
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FIG. 8: The scatter plot of ir + n~ invariant mass versus the 
Kg decay length for ip(2S) — > , yyK s K + iT~ + c.c. candidate 
events. 



After requiring \m n + n - — m K o\ < 0.015 GeV/c 2 , 
the 77 invariant mass is shown in Fig. 9 (a), and a 
clear ir° signal is seen. After requiring |m 77 — m 7r °\ < 
0.03 GeV/c 2 , the 7r ± 7r° invariant mass is shown in Fig. 9 
(b), where there is a clear p^ signal. 

The 77 invariant mass distribution is fitted with a 7r° 
signal shape determined with MC simulation plus a sec- 
ond order polynomial for the background, and the result 
is shown in Fig. 9 (a). The number of 7r° signal events 
fitted is 361 ±25, and the efficiency determined from MC 
simulation is 4.40%, including the effect of the interme- 
diate K ± p*K ( j^ state. 

The 7r^7r^ invariant mass distribution is fitted with a 
p ± signal shape determined with MC simulation plus a 
second order polynomial for the background, and the fit 
result is shown in Fig. 9 (b). The number of p ± signal 
events is 100±20, and the detection efficiency is 3.80% de- 
termined from MC simulation. The branching fractions 
of these two processes are determined to be 



B(ip{2S) -» 7r°K s K ± tt^) = (8.9 ± 0.6) x 10~ 4 , 
B{ip{2S) -» K ± p T K° s ) = (2.9 ± 0.6) x 10~ 4 , 

where the errors are statistical. 



FIG. 9: Invariant mass spectra of (a) 77 and (b) n ± Tv" for 
ip(2S) — > 77_ft'§.ff ± 7r T candidate events. Dots with error 
bars are data, the histograms are the fits using signal shapes 
determined from Monte Carlo simulation and second order 
polynomials for background, and the dashed curves are the 
background shapes from the fit. 



2. Signal analysis 

The event selection is similar to tp(2S) —> 
Tr°K s K + Tr~ +c.c, but only one photon is required. After 
event selection, the 7r + 7r~ invariant mass distribution is 
shown in Fig. 10. A fit is performed with a histogram de- 
scribing the K s shape obtained from MC simulation, the 
normalized histogram for ip{2S) — > ir a K s K + ir~ + c.c. 
background, and a Legendre polynomial for the other 
smooth backgrounds. The fit yields 115 ± 16 events. The 
detection efficiency is 4.83%, and the branching fraction 
is determined to be: 



B(tp(2S) -» "/KgK+Tr- + c.c.) = (25.6 ± 3.6) x 10~ 5 , 

where the error is statistical. Figure 11 shows the 
K%K 
tion. 



K s K ±r KT invariant mass distribution after event selec 



D. tp(2S) -f -yK+K-TT+n- 

For 4>{2S) — > -fK + K ~n + ir~ , candidate events require 
four charged tracks, among which two tracks must be 
identified as pions and the other two tracks identified as 
kaons. The background from ip(2S) — > Tr + Tr~J/ip is re- 
jected by requiring \m*e£,u ~ m JI^>\ > 0-05 GeV/c 2 , the 
background from "4>(2S) — > 72(7r + 7r _ ) is rejected by re- 
quiring X^K+K-n+n- < X 7 2( 7r + 7r -)> and thc background 
from tp(2S) — > -fK s K + n~ + c.c. is rejected by requiring 
\m w +Tc- - m K o \ > 0.04 GeV/c 2 . 

Figure 12 shows the K + K~Tr + iT~ invariant mass dis- 
tribution after event selection, where 361 events are ob- 
served. The backgrounds mainly come from tp(2S) — » 
TT°K + K~Tr + Tr~ final states including intermediate states. 
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FIG. 10: The -k + tt~ invariant mass distribution for ip(2S) — > 
'yKgK + ir^ + c.c. candidate events. Dots with error bars are 
data. The histogram is the fit with a histogram describing the 
Kg shape obtained from MC simulation, the normalized his- 
togram for ip(2S) — > tt°K s K + it~ + c.c. background (dashed 
histogram), and a Legendre polynomial for the other smooth 
backgrounds (dotted histogram). 
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FIG. 11: The KgK tr^ invariant mass distribution for 
ip(2S) — > ~(KgK + ir~ + c.c. candidates (dots with error bars). 
The shaded histogram is the sum of backgrounds including 
ip(2S) — > -k°KsK + tt~ + c.c. and continuum background. 



Using the x 2 fitting method of Section IV, the number 
of signal events is 132±19. The detection efficiency deter- 
mined from MC simulation is 4.94%, and the branching 
fraction for this process is determined to be: 



B(ip(2S) -» 1 K + K-Ti + Ti-) 
where the error is statistical. 



(19.1 ±2.7) x 10" 



E. i>(2S) -f -yK*°K+tt- + c.c. and 
4>(2S) -> -yK"°K*° + c.c. 

We apply the same event selection criteria as for 
ip(2S) — ► -/K + K ~ir + ir~ . For this decay channel, 
the main background channels are from ip(2S) — > 
K*°K-tt+tt° + c.c. (phase space), ip(2S) -» K* R-p + + 
c.c, and V(25) -» K*°K*°tt°. 



FIG. 12: The K + K~-k + -k~ invariant mass distribution for 
ip{2S) — > K~tt + it~ candidates (dots with error bars). 
The shaded histogram is background which mainly comes 
from iP(2S) -> it K+ K~ -K+ -k~ . 



Using a Brcit-Wigncr function to describe the signal 
along with the sum of normalized histograms from the 
V>(2S) -» K*°K-tt + tt° + c.c. (phase space), i>(2S) -» 



K*°K-p+ + c.c, and i/>(2S) 



K* K*°n background 
channels, and a second order Legendre polynomial to de- 
scribe other remaining backgrounds to fit the K^ir^ in- 
variant mass spectrum, 237±39 ip(2S) -» jK* K~-k + + 
c.c. candidate events are obtained, as shown in Fig. 13. 
Events from the intermediate state ip{2S) — ► ^K*°K* 
are counted twice, so an efficiency correction must be 
made for this. 
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FIG. 13: The K+tt* invariant mass distribution for tp(2S) — > 
■yK + K~tt + it~ candidates. Dots with error bars are data. 
The blank histogram is the result of a fit using a Breit-Wigner 
function to describe the signal along with the sum of normal- 
ized histograms from the tp(2S) — > K*° K~ -k + -k° +c.c. (phase 
space), ip(2S) -> K*° K~ p + + c.c., and ip(2S) -> K*°K*°-k° 
background channels, and a second order Legendre polyno- 
mial to describe other remaining backgrounds. The dashed 
histogram is the fitted sum of all backgrounds. 

The scatter plot of K + ir~ versus K~tt + invariant mass 
is shown in Fig. 14, where clear K*° and K signals are 
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seen. The numbers of ip(2S) — > ^K* Q K events and 
background events are estimated from the scatter plot. 
The signal region is shown as a square box at (0.896, 
0.896) GeV/c 2 with a width of 60 MeV/c 2 . Backgrounds 
are estimated from sideband boxes, which are taken 60 
MeV/c 2 away from the signal box. Background in the 
horizontal or vertical sideband boxes is twice that in the 
signal region. If we subtract half the number of events 
in the horizontal and vertical sideband boxes, we double 
count the phase space background. Therefore the back- 
ground is one-half the number of events in the horizontal 
and vertical boxes plus one fourth the number of the 
events in the four corner boxes. After subtraction, 41 ± 8 

ip(2S) -» ~/K*°K candidates are obtained, and the ef- 
ficiency is (2.75 ± 0.06)%. In simulating signal channels 
containing K* , the shape of K* is described by a P-wave 
relativistic Breit-Wigner, with a width 



ra 1 



r 2 p 2 



m 1 + 



y»2 ^2 



where m is the mass of the Kir system, p is the momen- 
tum of the kaon in the Kir system, T is the width of the 
resonance, m is the mass of the resonance, po is the mo- 
mentum evaluated at the resonance mass, r is the interac- 

l+r 2 2 

tion radius, and represents the contribution of the 

barrier factor. The value r = (3.4±0.6±0.3)( GcV/cT 1 
measured by the K~ir + scattering experiment [16] is used 
as an approximate estimation of the interaction radius r. 
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FIG. 14: The scatter plot of K + tv versus K 7r + invariant 
mass of ip(2S) — > ~{K + K~ -k + tt~ candidate events. The center 
box indicates the signal region for tP(2S) -> ~/K*°K events, 
and the other boxes are used for background determination. 

Taking into consideration the effect of the intermediate 
channel, i>(2S) -> -fK*°K*°, the efficiency for V(2S) -» 



jK*K 7T+ + c.c. is 6.86%, and we obtain the branching 
fractions: 



B(ip{2S) -» 1 K*K-tt+ + c.c.) = (37.0 ± 6.1) x 10" 5 
Bty{2S) -» ~fK*°K*°) = (24.0 ± 4.5) x 10~ 5 , 



where the errors arc statistical. 



F. ip(2S) ->■ j2(K+K~) 

For ip(2S) — > 'j2(K + K^), the candidate events must 
have four charged tracks, and every track must be iden- 
tified as a kaon. The backgrounds from tp(2S) — » 
72(7r + 7r~), jir + ir~ K + K~ are rejected by requiring the 
xic f° r ^ ne signal channel to be less than those for back- 
grounds. There are 15 events observed after event se- 
lection, and the 2(K + K~) invariant mass distribution is 
shown in Fig. 15. The detection efficiency for this channel 
is 2.93%. 
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FIG. 15: Invariant mass distribution of 2{K + K~) for 
, <p(2S) — > r y2(K + K~) candidates (dots with error bars). 
The shaded histograms is background mainly from tp{2S) — > 
tt°2(K+K-). 



The dominant background comes from ip{2S) — > 
w°2(K + K~). Using the branching fraction measured by 
CLEO [17], the estimated number of background events 
remaining is 8 ± 2. To measure the continuum contri- 
bution in this channel, the continuum data at E cm = 
3.65 GeV is analyzed using the same criteria as for ip(2S) 
data, and no events survive. It is also found that no 
events survive from the simulated 14 million inclusive 
tp{2S) decay MC sample. The upper limit on the num- 
ber of i>(2S) -» j2(K+K-) events is 14 at the 90% C.L., 
and the corresponding upper limit on the branching frac- 
tion after considering systematic uncertainties is 

B(ip{2S) -> j2(K+K-)) < 4.0 x 10~ 5 . 
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G. ip(2S) — ► 77T + 7r pp 

For tp(2S) — > jir + ir~pp, there must be four good 
charged tracks, and two of them must be identified 
as a proton anti-proton pair. The backgrounds from 
72(7T + 7r~) and jir + ir~K + K~ are rejected by requir- 
ing xlc f° r t ne signal channel to be less than for the 
background channels. To eliminate possible contamina- 
tion from tp(2S) — > n + -K~ J/tp, J /if) — > 7pp, we require 
Ke + c^ - > 0.02 GeV/c 2 . 

Figure 16 shows the ir + ir~pp invariant mass distribu- 
tion with 55 events after event selection. The detection 
efficiency for this channel is 4.47%. 




satisfying \m1^ oil — mj/^\ < 0.05 GeV/c 2 or hav- 
ing a pion pair in the K s mass region from 0.47 to 
0.53 GeV/c 2 . The remaining backgrounds mainly come 
from processes with multi-photon final states, such as 
tp{2S) -» 7r°3(7r+7r-),77r°3(7r+7r-), and 7r°7r°3(7r+7r-). 
Their contaminations are estimated using MC simula- 
tion. Figure 17 shows the 3(7r + 7r~) invariant mass dis- 
tribution after event selection with 118 events observed. 
The detection efficiency for this channel is 1.97%. Us- 
ing the x 2 fitting method described in Section IV, the 
upper limit on the number of signal events is 45 at 
the 90 % C.L., and the upper limit on the branch- 
ing fraction after considering systematic uncertainties is 
B(tp(2S) -> 73(7r+7r-)) < 17 x 10~ 5 . 
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FIG. 16: Invariant mass distribution of -k + tt~pp for tp(2S) — > 
77r + 7r~pp candidates (dots with error bars). The shaded his- 
togram is background mainly from ip(2S) — > -k°tt + -k~ pp and 

tp(2S) — » 7T + 7T~ J/tp, J/tp — » 7Pp. 

The dominant backgrounds are tp(2S) — > ir a ir + ir~pp 
and background remaining from tp(2S) — > 
tt + tt^ J/tp, J/tp — > jpp. The detection efficiencies for 
these two background channels are determined by MC 
simulation to be 0.35% and 0.18%, respectively. For the 
first background channel, using B(tp(2S) — > 7r°7r + 7r _ pj5) 
measured by CLEO [17], the number of background 
events remaining is estimated to be 35.8 ± 3.5. Simi- 
larly, the estimated number of background events from 
tp{2S) -> ir+n-J/tp, J/tp -> 7pp is 1.7 ±0.4. Subtracting 
backgrounds, the number of tp{2S) — > "fir + Tr~pp events 
is 17 ± 7, and the corresponding branching fraction is 

B(tp(2S) -> 77r+7r-pp) = (2.8 ± 1.2) x 10~ 5 , 

where the error is statistical. 



H. V(25) -»■ 73(77+77") 

For tp(2S) — > 73(7r + 7r _ ), six charged tracks are re- 
quired. The backgrounds from tp(2S) — > -k + t:~J/i)) 
and ip(2S) — + i^gifgTr+Tr - are removed by eliminat- 
ing events having the recoil mass of any pion pair 



FIG. 17: The 3(7r + 7r~) invariant mass distribution for 
tp(2S) — > 73(7r + 7r~) candidates (dots with error bars). 
The shaded histogram is background mainly from pro- 
cesses with multi-photon final states, such as tp(2S) — + 
7r 3(7r + 7O,77r (, 3(7r+7O, and TrV^Tr+O- 



I. tP(2S) -> ■y2(Tr + n-)K + K- 

1. Background estimation 

First, background from tp(2S) —* rjj/tp, r\ — > 
77T+7T - , J/tp — > tt + it~K + K~ is rejected by requiring 
m 2 ( 7r + 7r -)/f+if- < 2.9 GeV/c 2 . The dominant back- 
grounds remaining are tp(2S) — > 7r°2(7r + 7r _ )i ; C + i ; ir _ and 
77r°2(7r + 7r~)i ; r + . Branching fractions for these are 
not currently available, so they are measured using our 
tp(2S) data sample. 

For tp(2S) -> 7r 2(7r+7r-).?r+Ff-, the number of good 
photons is required to be iV 7 = 2 or 3. A kinematic fit 
is performed under the tp(2S) — > 772(7r + 7r _ )X + X _ hy- 
pothesis running over all selected photons, and the com- 
bination with the smallest \ 2 is retained. Background 
from tp{2S) -» ir+n-J/tp, J/tp -» 772(7r+7r-)if 
is rejected by requiring |m^ eC g i; — mj/^\ > 0.05 
GeV/c 2 . The possible backgrounds from tp(2S) — > 
72(7r + 7r _ )i ; 4' + i ; C _ and 372(7r + 7r~)i ; r + JT _ are rejected by 
requiring the x 2 value for the signal to be less than 
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those for the backgrounds. To remove backgrounds from 
4>(2S) — > jXcj decays, we require m 2 ^ +7t - )K+K - < 3.38 
GcV/c 2 . 

Eight main peaking background channels, including 
ip(2S) — > r)J/ip, TT°J/ip and 7Xcj to decay into the same 
final states, are simulated and fitted using the same pro- 
cedure and selection criteria, and 9.8 ± 4.5 background 
events are obtained. Using the 14 million inclusive MC 
sample, 14.7 ±5.7 background events arc found, which is 
consistent with the simulation result within the statisti- 
cal error. 

Figure 18 shows the 77 invariant mass distribution for 
ip(2S) -> tt°2(tt + tt-)K+K~ candidate events. A fit is 
performed with a n° signal shape determined from MC 
simulation plus a third order polynomial for the back- 
ground, and the number of ir° signal events is determined 
to be 57.4 ±9.8. After subtracting peaking backgrounds, 
the number of signal events is 47.6 ± 10.8. 

The detection efficiency determined from MC simu- 
lation using a phase space generator is 0.41% below 
m 2(-K+Tr-)K+K- = 3.38 GeV/c 2 , and the branching frac- 
tion is determined to be: 

B(tp(2S) -» 7r°2(7r + 7r- )K + K~) = (8.39 ± 1.91) x 10~ 4 

with the requirement m 2 ( 7T + 7r -)K+ k- < 3.38 GeV/c 2 , 
where the error is statistical. The effect of possible inter- 
mediate resonances is not considered. Assuming phase 
space production, the branching fraction extrapolated to 
the full m 2 r w + w -\ K + K - energy region is determined to be 
B{ip{2S) -» 7r°2(7r + 7r-)^+^-) = (11.5 ± 2.6) x lO" 4 . 
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FIG. 18: The 77 invariant mass distribution for tp(2S) — > 
'yy2(-ir + -K~)K + K~ candidate events. Dots with error bars are 
data, and the blank histogram is the fit with a 7r° signal shape 
determined from MC simulation plus a third order polynomial 
for the background. The curve is the fitted background. 

For i>(2S) -» 77r°2(7r+ n~)K + K-, the number of 
good photons is required to be JV 7 = 3 or 4. 



A kinematic fit is performed under the ip(2S) — > 
3"{2(Tr + Tr~)K + K~ hypothesis running over the se- 
lected photons; the combination with the smallest x 2 
is retained. Possible backgrounds from ip(2S) — ► 
(n'y)2(Tr + ir~)K + K~ with n = 1 2 and 4 and 
from ij){2S) -> 373(7r+7r-), 3jtt+tt-2(K+K-) 7 and 
3^K ± tt t 2(tt + tt~) are rejected by requiring that the 
X 2 of the signal is less than those of the back- 
grounds. Background from ip(2S) — > tt + tt~ J ftp, J/tp — > 
3jtt+tt~ K + K~ is rejected by requiring |TOrecoiZ — 
m J/-4>\ > 0-05 GeV/c 2 , and backgrounds from "4>(2S) — > 
7r°J/V>, rjJ/ip, ir°ir°J/iP, J/iP^2(tt+tt-)K+ K~ are re- 
jected with the requirement \m 2 ^ 7T + 7T -- )K + K - — m J/ip\ > 
0.05 GeV/c 2 . We select the 7r° from the 77 combi- 
nations as the one with m 77 invariant mass closest to 
m T o. To remove backgrounds from ip(2S) — > JXcj de- 
cays, m Tr o 2 ( Tr +Tr- )k+k- < 3.38 GeV/c 2 is required. 

After event selection, no significant 7r° candidates are 
observed. A fit with a n° shape determined from MC 
simulation plus a second order Legendre polynomial for 
background yields 27.1 ± 8.5 events. Fitting in the 
same way a histogram of 20 MC simulated background 
modes, 21.2 ± 10.6 background events are obtained. The 
detection efficiency determined from MC simulation is 
5.8 x 10 . The upper limit on the branching fraction at 
the 90% C.L., determined using POLE [18] and includ- 
ing systematic uncertainties, is 

B(ip(2S) -» jtt 2{tt + 7t-)K + K-) < 3.1 x 10~ 3 . 



2. Signal analysis 

For ip(2S) — > r )2(ir + TT~)K + K~ , six charged tracks are 
required, and two of them must be identified as kaons. 
The background from ip(2S) — > 7r + 7r~ J/ip, J/ip — > 
7± four charged particles is removed by requiring 
\ m recZii ~ m JlA > °-05 GeV/c 2 , and the back- 
grounds from ip(2S) — > 73(7r + 7r~), 7if ± 7T =F 2(7r + 7r~) and 
^2{K + K~ it~ are rejected by requiring the x 2 values 
for the signal to be smaller than for the backgrounds. The 
background from ip(2S) — > i)J/ip — > j2(n + ir~)K + K ~ is 
rejected by requiring m 2 ^+K-)K+K- < 2.9GeV/c 2 . 

Figure 19 shows the 2(tt + tt^)K + K ~ invariant mass 
distribution, where the shaded histogram is back- 
ground mainly from ip(2S) — » ir°2(ir + ir~)K + K~ 
and 77r°2(7r + 7r-)^+^-. For the ip(2S) -> 
n°2{-K + iT~)K + K~ background channel, the branch- 
ing fraction of (8.39 ± 1.91) x 10~ 4 is used since the 
MC sample is produced with rn 2 (. n +-K-)K+K- < 3.38 
GeV/c 2 . After subtracting all backgrounds, 17 events 
are obtained. The detection efficiency for this channel 
is 0.69%. The upper limit on the number of signal 
events is 15.5 at the 90% C.L., and the branching 
fraction after considering systematic uncertainties is 
B(ip(2S) -» 1 2{tt+tt-)K+K-) < 22 x lO" 5 . 
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TABLE I: Summary of systematic errors (%), where WR, e 7 , K% rec., and MC denote the wire resolution, photon efficiency, 
the error for K s reconstruction, and MC statistics, respectively. The sixth column gives the uncertainties due to the fits or 
the K* fit. 



Mode 


WR 


e 7 


PID 


Kg rec. 


fit 


Branching Fractions Background N^^s) 


MC Total 


1PP 


6.3 


2.0 


4.0 






— 


9.4 


4.0 


0.5 


12.8 




JppTT 


11.6 


6.0 


4.0 






— 


14.3 


4.0 


3.0 


20.4 


72(7T + 7T~) 


5.0 


2.0 


8.0 




3.0 (x fit) 


— 


6.4 


4.0 


1.0 


12.7 


■jKgK + TT + C.C. 


5.0 


2.0 




3.4 




— 


11.8 


4.0 


1.0 


14.1 


-fK + K 7T+7T 


10.7 


2.0 


8.0 




3-0 (x 2 fit) 


— 


17.2 


4.0 


2.2 


22.6 


r .j.n v *. __i_ _ 

7if °if + 7r + c.c. 


10.7 


2.0 


8.0 


— 


8.8 (K* fit) 


— 


10.0 


4.0 


1.2 


19.5 


-yK*°K*° 


10.7 


2.0 


8.0 


— 


— 


— 


10.0 


4.0 


1.1 


17.4 




10.4 


2.0 


8.0 


— 


— 


— 


20.6 


4.0 


1.1 


24.9 


72(7^+^-) 


11.1 


2.0 


8.0 


— 


— 


— 


14.2 


4.0 


1.6 


20.3 


73(7r + 7r _ ) 


5.0 


2.0 


— 


— 


— 


— 


— 


4.0 


2.2 


7.0 


72(^+71-)^+^" 


8.7 


2.0 


4.0 


— 


3.0 (x 2 fit) 




25.0 


4.0 


3.0 


27.5 




9.1 


6.0 


4.0 








6.0 


4.0 


3.1 


14.0 


pp-K^TT 


11.7 


8.0 


4.0 








4.3 


4.0 


1.7 


15.9 




10.0 


4.0 


8.0 








1.8 


4.0 


1.0 


14.2 


LU7V + 1T~ 


10.0 


4.0 


8.0 






0.8 (lu Br) 


2.0 


4.0 


1.0 


14.2 


w/ 2 (1270) 


10.0 


4.0 


8.0 






3.1 {h Br) 


10.0 


4.0 


1.0 


17.5 




10.0 


4.0 


8.0 






0.8 (6i Br) 


1.0 


4.0 


1.0 


14.1 




10.0 


4.0 




3.4 






3.0 


4.0 


1.5 


12.5 


K±p*K% 

TT° 2(7V+ TV' )K+K~ 


10.0 


4.0 




3.4 






8.0 


4.0 


1.5 


14.5 


13.5 


4.0 


4.0 








13.3 


4.0 


1.1 


20.2 




FIG. 19: The 2(ix + tt~)K + K~ invariant mass distribution for 
ip(2S) — > ~(2(ir + 7v-)K + K~ candidates (dots with error bars). 
The shaded histogram is background mainly from tp(2S) — > 
ir°2(-K+-K-)K+K- and j-K°2(n + -w-)K+ K~ . 

VI. SYSTEMATIC ERRORS 

Systematic errors on the branching fractions, listed in 
Table I, mainly originate from the MC statistics, the 
track error matrix, the kinematic fit, particle identifica- 
tion, the photon efficiency, the \ 2 fit method, the uncer- 
tainty of the branching fractions of intermediate states 
(taken from the PDG [3]), the uncertainty of the back- 
ground estimation, and the total number of ip(2S) events. 

1. The systematic error caused by the MDC tracking 
and the kinematic fit is estimated by using simula- 



tions with different MDC wire resolutions [10]. The 
systematic error ranges from 5% to 13.5% depend- 
ing on the number of charged tracks in the different 
channels. 

2. The photon detection efficiency was studied with 
J/ip — > 7T + 7r~7r° events [10], and the difference be- 
tween data and MC simulation is about 2% for each 
photon. 

3. Pure 7T and K samples were selected, and the parti- 
cle identification efficiency was measured as a func- 
tion of track momentum. On the average, a 1.3% 
efficiency difference per 7r track and a 1.0% differ- 
ence per K track are observed between data and 
MC simulation. We take 2.0% for each charged 
particle identification as a conservative estimate of 
the systematic error. 

4. In order to estimate the systematic error caused 
by the differences of the x 2 distributions between 
data and MC simulation, we use selected sam- 
ples of ip{2S) — > 7Xco, Xco — > K+K-tt+tt- and 
7T + 7r~7r + 7r~ to compare the \ 2 shapes of data and 
MC, because these two samples have similar final 
states and sufficient statistics. The difference is 
about 3%, which is taken as the systematic error 
of the x 2 fit method. We also performed an input- 
output study of the \ 2 fit, and found the differ- 
ence between input and output values is very small 
(< 0.5%) and is neglected. 

5. The background uncertainties are estimated by 
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changing the order of the polynomial or the fitting 
range used. The errors on the branching fractions 



of the main backgrounds (0(25) 



hs) have 



also been considered and included. The uncertainty 
of the background estimation varies from l%-25% 
depending on the channel and background level. 

6. The uncertainty of the total number of tp(2S) 
events is 4% [8] . 

Adding up all these sources in quadrature, the total 
systematic errors range from 7% to 28% depending on 
the channel. 



VII. RESULTS AND CONCLUSIONS 

Figure 20 shows the differential branching fractions for 
i/j(2S) decays into ~fpp, 72(7r + 7r~), r yK + K^7r + Tr~ , and 
r yKgK + ir^ + c.c, and the numbers of events extracted 
for each decay mode with mts < 2.9 GeV/c 2 are listed 
in Table II. Broad peaks, which are similar to those ob- 
served in J/ip decays into the same final states [14, 19], 
appear in the m P p and distributions at masses be- 
tween 1.9 and 2.5 GeV/c 2 and 1.4 and 2.2 GeV/c 2 , 
respectively. Possible structure within these broad peaks 
cannot be resolved with the current statistics. No ob- 
vious structure is observed in other final states. The 
branching fractions for < 2.9 GeV/c 2 in this paper 
sum up to 0.26% [20] of the total ip(2S) decay width, 
which is about a quarter of the total expected radiative 
-0(25*) decays. This indicates that a larger data sample 
is needed to search for more decay modes and to resolve 
the substructure of 0(25) radiative decays. 

Table III lists the results of '0(25') decays into 7r° + 
hadrons together with the world averaged values [3] , and 
values of Q h [= B{ip{2S) -> h)/B{J/^ -> h)]. For 
0(25) — > ir°2(ir + ir~) decay, intermediate resonances in- 
cluding a [/o(600)], / 2 (1270), w, and &i(1235) are ob- 
served, and the measurement of £>[0(25) — > cj/ 2 (1270)] 
agrees with the previous measurement using the same 
data sample [21]. The p ± resonance is observed in 
V>(25) -> n K^K + ir~ + c.c. decay mode. 

In summary, we report measurements of the branch- 
ing fractions of "0(25) decays into jpp, j2(jr + iT~), 
-fK^K+ir- + c.c, jK+R-tt+tt-, ^K*°K-t: + + c.c, 

-fK*°K*°, -/TT+n'pp, j2(K + K~), ^{tt+ti-), 
r )2{iv + -k~)K + K~ and the differential branching fractions 
for 0(25) decays into jpp, j2(tt + tt~), jK + K~tt + tt~ , 
and jKgK + n~ + c.c. with hadron invariant mass less 
than 2.9 GeV/c 2 . We also report branching fractions of 

_ n _ n 

ppn 7r 



■0(25) decays into 'yppir , pptt q tt°, it K%K + it -- <■■<■ 



K 



ir°2{n + iT-)K+K- 



The measurements of 0>(25) decays into tt°2(tt 1 r; 



and jTT 2{n + n-)K + K-. 

lott + tt , w/2(1270), and bfir T are consistent with 
previous measurements [3] and the recent measurements 
by the CLEO collaboration [17]. 



> 5 

CD 



O 

o 



o 

- 5 

£ 
— 
\ 

Q_ 



2.5 



TPP 



! 



yi\ TT 7T 7T 



++ + + 



77T + 7T K + K 



7KsK + 7T" + C.C. 



1 1.5 2 2.5 3 

m hs (GeV/c 2 ) 



FIG. 20: Differential branching fractions for 0(25) decays 
into 7pp, 72(7r + 7r - ), jK + K~ir + ir~ , and jKgK + -K~ + c.c. 
Here m^s is the invariant mass of the hadrons in each final 
state. For each point, the smaller vertical error is the statis- 
tical error, while the bigger one is the sum of statistical and 
systematic errors. 
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TABLE II: Results for tp(2S) — > 7 + hadrons. For each final state, the following quantities are given: the number of events in 
tp(2S) data, N Tot ; the number of background events from i/j(2S) decays and continuum, N Bg ; the number of signal events, N St9 ; 
and the weighted averaged efficiency, e; the branching fraction with statistical and systematic errors or the upper limit on the 
branching fraction at the 90% C.L. For all the radiative channels, except the jppiv and ■yiv 2(TT + iv~)K + K~ modes we require 
m h3 < 2.9 GeV/c 2 . The branching fraction for yn°2(iv + -K~)K + K~ is measured with the requirement rn lr o 2 (i r +n-)K+K- < 3.38 
GeV/c 2 . Possible interference effects for the modes with intermediate states are ignored. 



Mode 


jyl ot 


N B9 




e(%) 


Z3(xl0 b ) 


1VP 


ozy 


1 07 

lo ( 


1 A Q _l_ 1 O 


OO.O 


o_i_n /i_i_n a 
z.y±U.4±U.4 


7PP7T° 


345 


219 


126 ±38 


8.94 


10.1 ±3.1 ±2.1 


72(7T- + 7T~) 


1697 


1114 


583 ±41 


10.4 


39.6±2.8±5.0 


-yK^K+TT- ± c.c. 






115 ± 16 


4.83 


25.6±3.6±3.6 


■yK+R-n+TY- 


361 


229 


132 ± 19 


4.94 


19.1±2.7±4.3 


■yK*"K + n- + c.c. 






237 ± 39 


6.86 


37.0±6.1±7.2 


jK*°K*° 


58 


17 


41 ±8 


2.75 


24.0±4.5±5.0 




55 


38 


17 ±7 


4.47 


2.8±1.2±0.7 




15 


8 


< 14 


2.93 


< 4.0 


73(7r + 7r _ ) 


118 


95 


< 45 


1.97 


< 17 


j2(tt+tt-)K+K- 


17 


13 


< 15.5 


0.69 


< 22 


jtt"2(tv + tt-)K+K- 


27 


21 


< 24.9 


0.058 


< 310 



TABLE III: Results of ip(2S) — > tv° ± hadrons. Here N St9 is the number of signal events, e is the detection efficiency, B is 
the measured branching fraction, £! PDG is the world averaged value [3], and Qh = B(ip(2S) —* h)/B(J/tp — > h). The branching 
fraction for 7r°2(7r + ~tt~)K + K~ is measured with the requirement m,2( n + w -)K+K- < 3.38 GeV/c 2 . 



Mode: h N»' 9 g(%) g(x!Q- 4 ) ^""'(x 10~ 4 ) Q h (%) 

7r u 2(7r+7r-) 2173 ±53 6.32 24.9 ± 0.7 ± 3.6 23.7 ± 2.6 10.5 ± 2.0 

ujtv+tv- 386 ± 23 3.74 8.4 ±0.5 ±1.2 6.6 ±1.7 11.7 ±2 .4 

w/ 2 (1270) 57 ± 13 3.65 2.3 ± 0.5 ± 0.4 2.0 ± 0.6 5.4 ± 0.6 

bfiv^ 202 ± 21 3.24 5.1 ± 0.6 ± 0.8 3.6 ± 0.6 17.0 ± 4.2 

PPty°ty° 203 ± 27 8.30 1.75 ± 0.21 ± 0.28 

■k ^^^ 361 ± 25 4.40 8.9 ±0.6 ±1.1 

if ± p =F Kg 100 ± 20 3.80 2.9 ±0.6 ±0.4 

ir°2(ir + ir-)K+K- 48 ±11 0.41 8.4 ±1.9 ±1.7 
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